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М?’яка математична морфологя 
для фалльтраци б1нарних зображень 


В статье рассмотрены операторы мягкой морфологии и алгоритмы удаления шума с бинарных изображений 
на их основе. Алгоритмы отличаются наличием таких дополнительных параметров работы, как порог 
фильтрации и количество итераций, а также возможностью удаления шума вида островки, проколы, 
выступы и вырывы за один проход. Применение данных алгоритмов позволяет привнести меньше 
искажений при работе, даже если объекты изображения располагаются на расстоянии 1 пикселя друг 
от друга. Приведены результаты экспериментальной проверки работы на примере изображения 
фотошаблона печатной платы. 
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Введение 


Базовые операторы классической математической морфологии — эрозия и дила- 
тация [1-3] выражаются как 


Рисие(А,В) = 9 {А+Ь:Ь ЕВ}, (1) 


Ето4е(А,В) = п ЗА-Б:БЕВ}, (2) 
где А — исходное изображение, В — структурирующий элемент. 
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Для удаления шума с изображения могут применяться фильтры, основанные 
на операциях отмыкания и замыкания, которые представляют собой последовательность 
эрозии и дилатации 


Ореп(А,В) = ОИже(Ето4е(А,В),В). (3) 
С]о5е(А,В) = Егодае(ОИае(А,-В),-В). (4) 
В [4], [5] были предложены операторы мягкой морфологии, а также фильтры 


на их основе, — мягкая эрозия и дилатация сигнала /: 2" > по структурирующей си- 
стеме [В‚,В>,К] определяется как: 


те [В,В»2,К]| (х) = К-у наименьшему значению мультимножества 


(5) 
{КО ДЬ>) : БЕ В2(%) 9 {КБ : БЕ (Вх) для всехх Ее #", 


ТФ [ВьВ2,К | (х) = К-у наибольшему значению мультимножества 
{КОДЬ>2) : БЕ В>х)} 9 К) : В, Е (ВК) для всеххЕ 2". 
Структурирующая система [В‚,В›,^] состоит из 3 параметров: конечных мно- 


(6) 


жеств Вг и В›, В› Е В, при этом В‚,В› Е Де. натуральное число К удовлетворяет 


условию 1 < К < Саи(В/). Множество В = В! Ц В», В, П В? = © называется структу- 
рирующим множеством, В› — центром, В; — границей, А — порядковым индексом 
центра (параметр повторяемости). 

В [6], [7] предложено следующее определение операторов бинарной мягкой 
морфологии: 


Ае [В,В>А (>) = {хЕА| (Кх Са ЙА П (ВР + Са А П (В) (1) 
> кх Сага[В/] + Са В — +1, 


АФ [В,,В> К] (<) = {хЕА| (Кх СанйА П (ВР) + Са А П (ВХ > К, (8) 
где К — порядковый индекс, который определяет, сколько раз элементы ядра 


учитываются в конечном результате. Если КА = 1 или В = В, (В> = ©), то операторы 
мягкой морфологии превращаются в стандартные операторы математической мор- 
фологии. 

Операторы математической морфологии, основанные на использовании нечет- 
ких множеств, предложены в [8]. В данном подходе нечеткость определяется тем, 
насколько структурирующий элемент вписывается в изображение. На сегодняшний 
день не существует единого подхода для определения нечеткой эрозии и дилатации. 
В [9-18], [8] даются различные определения базовых операторов нечеткой мор- 
фологии, отличающиеся результатами работы. Например, согласно [8], операторы 
нечеткой эрозии и дилатации по нечеткому структурирующему элементу в терминах 
функций принадлежности имеют следующий вид: 


Илев(х) = п (па [1, 1 + мл(х + У) - иво) 


9 
= ши], пуп [1 + мл(х + У) - ивОТ, йе 
Илев (х) = тах [пи [0, дл(х - у) + ивО>) - 1] 
УЕВ (10) 


= тах[ 0, тах [ма(х - у) + ив» - |, 


гдех, у Е 27 - пространственные координаты, а ил и ив— функции принадлеж- 
ности изображению и структурирующему элементу. 
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В [19-21] приводится описание подхода, объединяющего нечеткую и мягкую 
морфологию. Операторы мягкой нечеткой эрозии и дилатации выглядят следующим 
образом (определение для нечеткой морфологии использовалось согласно [8]): 


дл абы, вз, (%) = шит, тип, СКО (аб + У) - ив + 0} 5 


{ ма(х- 2) — ив2(2) + 1})], 


= @® ы — 
паев, во» (х) = тах[0, — шах (КО (их) тив - 0} 9 (12) 


{ ма(х—2) + ив2(2) — 1})], 


о = а пространственные координаты, а ил, ив/ и ив? — функции принад- 
лежности изображению А, ядру В! и мягкой границе В? структурирующего элемента. 


(0 


Для нечеткого структурирующего элемента выполняются условия В Е 2 АВ 


В, П В2 = ©. Если К = 1, то операторы мягкой нечеткой морфологии превращаются 
в операторы нечеткой математической морфологии. 

Следует отметить, что публикации в литературе подтверждают актуальность 
использования мягкой [22-25] и нечеткой [15], [26-28] морфологии для обработки изобра- 
жений. 


Алгоритмы фильтрации изображения на базе 
операторов мягкой морфологии 


В [29], [30] предложены операторы мягкой морфологии, которые сочетают в 
себе элементы мягкой и нечеткой морфологии: 

1) как и в мягкой математической морфологии, предлагаемый структурирующий 
элемент состоит из мягкой границы, но отличается тем, что ядро структурирующего 
элемента равно пустому подмножеству; 

2) как и в нечеткой морфологии, принцип работы основан на возможности вписать 
структурирующий элемент в изображение, но сам элемент не рассматривается как не- 
четкое множество. 

Термин «мягкий» используется для того, чтобы показать, что для одного и того 
же структурирующего элемента можно получить различные результаты, изменяя 
значение порога фильтрации. 

Операторы мягкой морфологии имеют следующий вид: 


п.+1<п, >а=0 

5оПЕто4е( А, В!) =4 ^ ‚аЕА,, (13) 
п, +Ё> п, >а=1 
п, +1<и, >а=1 

бойрИае( А, В,Г) =3 ‘ ‚ЕЛ, (14) 
п, +Ё> т, >а=0 


элемента, совпавших по значению с ненулевыми пикселями структурирующего 
элемента, и, =У(а, лЬ,) — число пикселей на изображении в маске структурирующе- 
7 


го элемента, значение которых не совпало со значением соответствующих пикселей 
структурирующего элемента, + — порог фильтрации, а — пиксели исходного черно-белого 
изображения 4, 6 пиксели плоского структурирующего элемента В. Термин — плоский 
структурирующий элемент обозначает, что его пиксели могут принимать значения 0 и 1. 
Исходными данными является бинарное изображение. Параметры работы — 
размер и форма структурирующего элемента, значение порога фильтрации или зна- 
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чение порогов для операций эрозии и дилатации, если для этих операций используются 
разные пороги. Результат — бинарное изображение. 

Для удаления шума с изображения предлагается использовать операции мягкого 
отмыкания и замыкания, которые имеют вид: 


5о{Ореп(А,В,1) = ЗойОИеие(5ойЕто4е(А, В,1),В,1), (15) 


5оС1озе(А, В,1) = ЗойЕто4е(ЗойОИеие(А,В,1),В,1). (16) 


Исходными данными является бинарное изображение. Параметры фильтрации — 
размер и форма структурирующего элемента, количество проходов фильтра, значение 
порога фильтрации или значение порогов для операций эрозии и дилатации, если 
для этих операций используются разные пороги. Результат — бинарное изображение. 

Мягкое отмыкание и замыкание не обладают свойством идемпотентности, поэтому 
могуг многократно применяться с одним структурирующим элементом, и это свойство 
предлагается использовать для конструирования мягких морфологических фильтров. 
Алгоритм фильтрации изображения на основе операций мягкого отмыкания и замы- 
кания приведен на рис. 1. 

Параметры фильтрации -— это структурирующий элемент В (его размер и форма), 
порог фильтрации { и количество итераций фильтра /. 

Также можно использовать последовательные мягкие морфологические фильтры, 
которые обозначаются 5ОС, 5СО, 5ОСО, 5СОС и т.д., где сокращениями 50 и 5С обозна- 
чаются операции мягкого отмыкания и замыкания. Данные фильтры сглаживают неров- 
ности краев элементов аналогично фильтрам классической математической морфологии. 
А с предложенными операторами мягкой морфологии можно использовать фильтры 
оз О, е5С, БобОБе5С, уе СТоЭО и Т.Д., ГДЕ 1» — количество итераций мягкого отмыка- 
НИЯ, а 1. — количество итераций мягкого замыкания. Количество проходов фильтра 
задается в зависимости от качества входного изображения. 


Начало 


1. А= бойрйаю(А.В) 1 -- ЗоНВтоЧеА.В) 


Параметры: 


Параметры: 
1 — Във 


Въ п 


2. А= ЗоЙЕюас(АВ) 2. А = бой ОЦа (А.В) 


3. Все итерации 


выполнены? 
а) 6) 


Рисунок 1 — Алгоритмы фильтрации операциями мягкой морфологии: 
а) отмыкание, 6) замыкание 


3. Все итерации 
выполнены? 
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Экспериментальные результаты применения 
мягкого морфологического фильтра 


На примере изображений фотошаблона ПП рассмотрим преимущества и особен- 
ности применения мягкого морфологического фильтра. На рис. 2а приведен фрагмент 
исходного полутонового изображения фотошаблона ПП с шумом. Размер изображе- 
ния 1596 х 1007 пикселей ‚ | пиксель изображения соответствует 42,3 мкм. 
Результат операции мягкого отмыкания по квадратному структурирующему элементу 
размером 3 х 3 пикселя с порогом 0 приведен на рис. 26, а результаты классических 
операций отмыкания и замыкания по тому же структурирующему элементу на рис. 
2в и г, соответственно. Оценка качества фильтрации вычислялась как 


О = 1 -— соци: (Тез ХОВ Ею/оп) / чате (Тез), (17) 
где Тез! — результат фильтрации, Е1оп — изображение эталонного фотошабло- 


на без шума, соип ) — операция вычисления количества ненулевых пикселей. Оценка 
приведена в табл. | и 2. 


в) 
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Рисунок 2 — Результаты работы морфологического фильтра по квадратному 
структурирующему элементу размером 3 х 3 пикселя: а) исходное изображение 
фрагмента фотошаблона печатной платы с шумом, 6) мягкое отмыкание 
ЗоЯОреп(А,В,0), в) отмыкание Ореп(А,В), г) замыкания С]озе (А,В), 

д) СозеОреп(А,В), е) ОрепСТюзе (А,В) 

Таблица 1 — Качество фильтрации при использовании структурирующего элемента 
размером 3 х 3 пикселя 


Значение порога Качество фильтрации 
фильтрации 1 ЗойОреп(А, В, 1) ЗойСо5е(А, В,1) 
В - квадрат В-ромб В - квадрат В-ромб 

0 0.99286 0.98481 0.99286 0.98876 
2 0.99172 0.9901 0.99116 0.99135 
4 0.99254 0.98795 0.99045 0.98795 
6 0.99167 0.99135 0.98832 0.9901 
8 0.99061 0.98876 0.98532 0.98481 

Классическая Ореп(А,В) Ореп(А,В) С1озе(А, В) Созе(А, В) 

морфология 0.99061 0.98876 0.98532 0.98481 


Таблица 2 — Качество фильтрации при использовании структурирующего элемента 
размером 5 Хх 5 пикселей 


Значение порога Качество фильтрации 
фильтрации # ЗойОреп(А, В,Е) ЗоЯС!озе(А, В,г) 
В - квадрат В-ромб В - квадрат В-ромб 
0 0.99006 0.98189 0.98972 0.98866 
2 0.98945 0.9889 0.98872 0.99209 
4 0.9885 0.98956 0.98756 0.99218 
6 0.98896 0.99011 0.98731 0.99195 
8 0.98998 0.99125 0.98771 0.99212 
10 0.99065 0.99263 0.988 0.99284 
12 0.9897 0.99428 0.98659 0.99428 
14 0.98855 0.99284 0.9854 0.99263 
16 0.98898 0.99178 0.98455 0.99096 
18 0.98928 0.99195 0.98417 0.99057 
20 0.98998 0.99218 0.98297 0.98993 
22 0.98866 0.99209 0.97961 0.98808 
24 0.98599 0.98866 0.97519 0.98157 
Классическая Ореп(А,В) Ореп(А,В) Созе(А,В) С озе(А, В) 
морфология 0.98599 0.98866 0.97519 0.98157 
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Классическая операция отмыкания удаляет островки с изображения; сглаживает 
выступы; увеличивает вырывы; образует разрывы. А классическая операция замыка- 
ния удаляет проколы с изображения; сглаживает вырывы; увеличивает выступы; обра- 
зует короткие замыкания. 

Для высокого качества фильтрации необходимо удалить как островки, так и проко- 
лы, сгладить как выступы, так вырывы. Поэтому можно применить последователь- 
ный фильтр С!озеОреп(А,В) (рис. 29) или ОрепСозе(А,В) (рис. 2е). Оценка качества 
фильтрации приведена в табл. 3 и 4. 


Таблица 3 — Качество работы последовательного фильтра с использованием 
структурирующего элемента размером 3 х 3 пикселя 


Значение Качество фильтрации 
порога ЗойОрепСТозе(А, В, 1) ЗойСю5еОреп(А, В, 1) 
фильтрации 1 В - квадрат В-ромб В - квадрат В-ромб 
0 0.99462 0.99086 0.99462 0.99067 
2 0.99101 0.99352 0.99127 0.99338 
4 0.99164 0.98804 0.99191 0.98804 
6 0.99356 0.99338 0.99363 0.99352 
8 0.99382 0.99067 0.99397 0.99086 
Классическая ОрепСозе(А,В) ОрепС1озе(А,В) СозеОреп(А,В) СозеОреп(А,В) 
морфология 0.99382 0.99067 0.99397 0.99086 


Таблица 4 — Качество работы последовательного фильтра с использованием 
структурирующего элемента размером 5 х 5 пикселей 


Значение Качество фильтрации 

порога ЗойОрепСТозе(А, В, 1) ЗойСю5еОреп(А, В, 1) 
фильтрации { В квадрат В-ромб В квадрат В -ромб 
0 0.98695 0.99092 0.98697 0.99105 
2 0.98416 0.99388 0.9836 0.99397 
4 0.98241 0.99005 0.98216 0.99023 
6 0.98218 0.98934 0.98253 0.98959 
8 0.98345 0.98969 0.98425 0.98887 
10 0.9867 0.99113 0.98662 0.99117 
12 0.98417 0.99512 0.98326 0.99512 
14 0.98315 0.99117 0.98254 0.99113 
16 0.98389 0.98837 0.98371 0.98902 
18 0.98565 0.99006 0.98485 0.98972 
20 0.98697 0.99078 0.98675 0.99041 
22 0.98799 0.99318 0.98778 0.99335 
24 0.98687 0.99084 0.98665 0.99082 
Классическая ОрепСТозе(А,В) ОрепСТо$е(А,В) СозеОреп(А,В) СозеОреп(А,В) 
морфология 0.98687 0.99084 0.98665 0.99082 


Классический последовательный фильтр лучше удаляет шум, чем фильтр, исполь- 
зующий только одну операцию вида отмыкание или замыкание. Последовательный 
фильтр на базе операторов мягкой морфологии лучше удаляет шум, чем последова- 
тельный фильтр, использующий операции классической мягкой морфологии. Одна- 
ко из рис. 26 видно, что фильтр, использующий только операцию мягкого отмыкания 
или замыкания, удаляет как островки, так и проколы, сглаживает и выступы, и вырывы 
за один проход, поэтому целесообразным является выбор именно такого фильтра в случае 
необходимости высокой скорости работы, т.е. когда незначительно повышение качества 
фильтрации (например, с 0.99286 до 0.99462, как в случае с операциями ЗоЙОрек(А, В, 1) 
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и $0 С! озеОреп (А, В,1)) является меньшим преимуществом, чем сокращение вычисле- 
ний в процессе фильтрации в 2 раза. 

На рис. 3 показан фрагмент фотошаблона с выделенными областями обрыва и ко- 
роткого замыкания. 


с” = = 
Рисунок 3 — Фрагмент фотошаблона с выделенными областями 
разрыва и короткого замыкания 


На рис. 2 видно, что применение операций классической морфологии для удаления 
шума привело к смыканию обрыва и разрыву короткого замыкания на исходном за- 
шумленном изображении. 

Если отфильтрованное изображение подать на вход системы контроля дефектов 
топологии, то данные дефекты не будут найдены. Отсюда следует недостаток фильтра, 
использующего классические операции отмыкания, замыкания или их комбинацию, 
который заключается в недостаточной чувствительности такого фильтра, то есть обра- 
зовании смыканий или обрывов на изображении при близком расположении объектов. 

В данном случае, под близким следует понимать расположение объектов на рас- 
стоянии меньше радиуса структурирующего элемента. Например, минимальным сим- 
метричным структурирующим элементом, использование которого не сдвигает изобра- 
жение, является ромб размером 3 Хх 3 пикселя. 

Радиус такого элемента равен 2, и при расстоянии между объектами, равному 1 
пикселю, происходит искажение топологии, т.е. образуются короткие замыкания или 
обрывы. 

Использование операторов мягкой морфологии с дополнительным парамет- 
ром - порогом фильтрации, позволяет сохранить топологию печатной платы в процессе 
удаления шума, даже если объекты изображения располагаются на расстоянии 1 пик- 
селя друг от друга. 


Выводы 


На базе операторов мягкой математической морфологии предложены алгорит- 
мы фильтрации, отличающиеся наличием дополнительных параметров работы, таких 
как порог фильтрации и количество итераций. В отличие от классических фильтров 
на базе отмыкания и замыкания, фильтр обладает свойством удалять шум вида островки, 
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проколы, выступы и вырывы за один проход без применения последовательных филь- 
тров. Применение данных алгоритмов позволяет привнести меньше искажений при рабо- 


те, 


даже если объекты изображения располагаются на расстоянии 1 пикселя друг от друга. 
Экспериментальная проверка работы мягкого морфологического фильтра на при- 


мере изображения фотошаблона печатной платы показала лучшие результаты по срав- 
нению с классическими морфологическими фильтрами. 
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вой Майетайса! Могрйо1озу юг ЕШептиз о} Втагу [тазез 

Тре рарег ргорозез а зой тогрбо|оз1са| орегаюг$ ап аеогИНил$ 10 гетоуе по15е 
Пот Ыпагу ппасез. Тре Ееабигез ое а1еогИ ит - 1$ Ше ргезепсе оРа@4опа| рагатеегз 
зиср аз Ге гезВо14 ог Него, питльег оЁ НегаНоп$, ап4 {йе ау {ю гетоуе по1зе Югт 
бе 131ап4$, рт Боез, тоизебйез апа зригз ш опе раз$. АррИсаНоп ог 'Пезе а1еотИНил$ сап 
тис [е5$ Ч1$ютИоп дигте орегаНоп, еуеп № \е ппазе об]ес{5 аге зрасе4 опе рхе| ара. 
Треге аге гези $ оРе ехрегипета1 уемйсаНоп оЁ Фе ехатр!е оРе РСВ ай\о!К. 
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